v.i from ideal and viscous hydrodynamic simulations of nuclear collisions at the BNL 
Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC) 



O 

(N 



o 

=5 



(N 
> 

(N 
O 
(N 

O 
O 



Matthew Luzum, Clement Gombeaud, and Jean- Yves Ollitrault 
Institut de Physique Theorique, 
CEA, IPhT, F-91191 Gif-sur-Yvette, France 
CNRS, URA 2306, F-91191 Gif-sur-Yvette, France 
(Dated: June 1, 2010) 

We compute V4/(v2) 2 in ideal and viscous hydrodynamics. We investigate its sensitivity to details 
of the hydrodynamic model and compare the results to experimental data from RHIC. Whereas i>2 
has a significant sensitivity only to initial eccentricity and viscosity while being insensitive to freeze 
out temperature, we find that V4,/(v2) 2 is quite insensitive to initial eccentricity. On the other hand, 
it can still be sensitive to shear viscosity in addition to freeze out temperature, although viscous 
effects do not universally increase Vi/{v2) 2 as originally predicted. Consistent with data, we find 
no dependence on particle species. We also make a prediction for Vi/(v2) 2 in heavy ion collisions at 
the LHC. 



I. INTRODUCTION 

Much can be learned about the bulk properties of the 
medium created in relativistic heavy ion collisions by 
studying the azimuthal distribution of emitted particles. 
Typically, the dependence on the azimuthal angle <fi with 
respect to the collision plane is written as a Fourier series 
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In a collision between identical nuclei at midrapidity, the 
sine terms and all odd terms are negligible when aver- 
aged over many events. There has been much study of 
the elliptic flow coefficient, v 2 [2], which has been an 
important observable for determining, e.g., the viscos- 
ity of the medium 0-H. Less studied, however is the 
next harmonic u 4 [7, 8]. In this article we use viscous 
and ideal hydrodynamic simulations to determine what 
«4 can tell us about the heavy ion collision fireball cre- 
ated at the Relativistic Heavy Ion Collider (RHIC), and 
to make predictions for planned heavy ion collisions at 
the Large Hadron Collider (LHC). 

For details of the hydrodynamic model used, the reader 
is directed to Ref. [2|. As a brief summary, we solve 2+1 
dimensional conformal viscous hydrodynamic equations 
(with zero bulk viscosity) using both Glauber- and CGC- 
type initial conditions and a Cooper-Frye freeze out rou- 
tine followed by resonance feed down. All calculations 
use the same parameters as those that gave the best fit 
to RHIC data [2J (see Table HJ), except where specified. 



II. REVIEW 

First recall the simplest prediction for V4 @ that pre- 
dicts V4 — 0.5(^2) 2 for a boosted thermal freeze out dis- 
tribution. To summarize, in the low temperature limit, 
one can solve the zero-viscosity Cooper-Frye freeze out 



Init. cond. 


v/s 


Ti [GeV] 


T f [GeV] 


T [fm/c] 


Glauber 


10" 4 


0.340 


0.14 


1 


Glauber 


0.08 


0.333 


0.14 


1 


Glauber 


0.16 


0.327 


0.14 


1 


CGC 


1Q -4 


0.310 


0.14 


1 


CGC 


0.16 


0.299 


0.14 


1 



TABLE I. Parameters used for the viscous hydrodynamics 
simulations except as indicated in each relevant section. The 
bold faced lines give the best fit to V2, (pt), and multiplic- 
ity 1. 



formula with Boltzmann statistics 
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by performing a saddle point approximation. At large 
transverse momentum (i.e., where the minimum of p-u > 
to), the dominant part of the integral comes from the 
part of the freeze out surface where the fluid velocity is 
parallel to the momentum of the emitted particle, and is 
at a maximum. We can decompose this maximum fluid 
4-velocity magnitude it max at each angle in a Fourier 
series 

«max(^) = U(l + 2V 2 cos(20) + 2 V 4 cos(4<£) + ...). (3) 

Plugging this into Eq. ([2]) in a saddle point approximation 
and expanding to leading order in V2 and V4, one obtains 
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with rat — V 'pt + m 2 and v = U/y/l + U 2 . At a large 
enough p t , the first term in Eq. ([S]) dominates, and it is in 
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this sense that it can be said that V4 is largely generated 
by an intrinsic elliptic flow (i.e., V2) rather than a fourth 
order moment of the fluid flow (V4). As one relaxes the 
large pt limit, the ratio of 1*4 to (2*2) 2 should have cor- 
rections from 1/2 that depend on the particular collision, 
but behave as 1/pt (or equivalently 1/1*2). 

This derivation is the motivation for studying the 
quantity 1)4/(112) 2 rather than V4 alone, and is the quan- 
tity to be studied here. 

Equation ([5]) is quite general and does not depend on 
the details of the freeze out surface or the fluid flow that 
produces it. It should be noted, however, that several 
assumptions must be made in order to derive it: small 
temperature, large transverse momentum, small V2 and 
1)4, and zero viscosity. As these restrictions are relaxed, 
the behavior becomes more complicated and more sensi- 
tive to details. Therefore, in order to fully understand the 
final results of the most realistic simulations as well as to 
make contact with previous results from ideal hydrody- 
namics and transport simulations, it will be useful to gain 
insight by starting with the simple case of vanishing vis- 
cosity and low freeze out temperature, and then adding 
viscosity and a realistic freeze out temperature once the 
generic features of the simpler case are understood. 



III. REVIEW OF RHIC DATA 

W4 has been measured at RHIC by both the STAR (Tol — 
GJ] and PHENIX [ll-GJ] collaborations. They both find 
a significant deviation of i*4/(i>2) 2 from a constant 1/2. 
Rather, the reported value of the ratio is close to a con- 
stant value of 1 for a large range of transverse momentum 
(p t > 500 MeV) Ql], centrality (~10-70%) HE El, and 
for all identified particles that have been measured (pi- 
ons, kaons, and protons [13, EH). ^ i s l ar g er for more 
central events, while there is no clear consensus for small 
Pt and very peripheral events. The results from PHENIX 
tend to be somewhat smaller than from STAR, likely 
indicating different sensitivities to non-flow effects, al- 
though all results are still well above 0.5. 

Much of this discrepancy between prediction and ex- 
periment can be understood by noting that both V4 
and V2 are averaged over many events before the ra- 
tio is taken. Even if the ratio 174/(1*2) 2 were exactly 
1/2 in each event, any fluctuation event- by- event (both 
at a given centrality and within a centrality bin) will 
tend to increase the measured 1*4/(1*2 ) 2 - In fact, using 
a Glauber Monte Carlo model of eccentricity fluctua- 
tions can explain most of the difference, as demonstrated 
in Ref. [2(| ■ With this understanding, all calculations 
here — since they lack fluctuations — should be understood 
as relating to measured data by an overall (momentum 
independent) multiplicative factor that depends on the 
centrality bin chosen. 



IV. RHIC RESULTS 

We begin by noting that, since 1*4 tends to be much 
smaller than V2 , numerical errors can be of greater con- 
cern. We estimated this uncertainty in two ways. First, 
we looked at the results at b = 0, where V4 should be zero 
theoretically. However, the use of a square lattice as well 
as the presence of numerical/rounding errors results in a 
calculated V4 that can be as much as ~ 5% of the value 
calculated at, e.g., b = 6 fm (or a much smaller fraction 
of the value for more peripheral collisions.) Second, we 
reran a number of simulations with half the usual grid 
spacing (0.2 fm instead of 0.4 fm). Similarly, we found 
that V4/(v 2 ) 2 changed by less than 5%-10% as long as 
we only considered collisions for b > 5-6 fm. Therefore, 
we only report the results for b = 6 fm and larger. Reli- 
able results outside this range can be obtained with more 
computer time or a better hydrodynamics algorithm (as 
reported in Ref. plj). 

In addition to numerical uncertainties, there are sev- 
eral other concerns to keep in mind. Although the results 
are numerically stable above p t ~0.2-0.3 GeV, they are 
moderately sensitive to details (e.g., small changes in pa- 
rameters or the neglect of resonance feed down) up to 
~0.5-l GeV, although they seem quite robust at larger 
p t . Of course, one must also note that in reality the hy- 
drodynamic description should break down at some large 
value of pt and therefore one has less and less confidence 
in the results as pt increases. However, it is not clear a 
priori where exactly this break down should occur, and 
so it is still useful to calculate quantities even at large 
Pt- Only comparison to data can then give information 
about this break down of the hydrodynamic description. 

None of these concerns should change the trends and 
conclusions reported here. 

Figure Q] shows the results for pions with 77/s = 0.0001 
(corresponding to ideal hydrodynamics), and all param- 
eters as specified in Table Q] except a small freeze out 
temperature of 100 MeV. As found previously in ideal 
hydrodynamic simulations [H, EH, ^4/(1*2) 2 approaches 
a constant value of roughly 1/2 at asymptotically large 
values of pt, and increases as 1/pt for decreasing p t . In 
contrast to expectations and data, however, there is a 
definite dependence on impact parameter. Part — but not 
all — of this dependence as well as the deviation from 1/2 
can be understood by the fact that here t>2 is unphysically 
large, as we have turned off viscous effects. For example, 
with CGC initial conditions at b = 11 fm and pt = 3.5 
GeV, v 2 reaches 70%. i*4/(?-*2) 2 would then be expected 
to reach almost 0.62 instead of 0.50 at large pt if one 
assumes V4 of the fluid to be negligible (see Ref. [20]). 

Figure [T] also makes it apparent that, unlike 1*2 itself, 
the results for i*4/(t>2) 2 are quite insensitive to the initial 
eccentricity. Although there is a dependence on impact 
parameter, it is small, and there is little difference for 
Glauber and CGC initial conditions, despite having a 
very different initial eccentricity (and therefore very dif- 
ferent V2 for the same viscosity — see, e.g., Figs. 4 and 8 
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Glauber, r|/s = 0.0001, T f = 100 MeV 



CGC, r|/s = 0.0001, T f = 100 MeV 





FIG. 1. (Color online) Pion V4,/(v2) 2 from ideal hydrodynamics calculations with Glauber and CGC initial conditions and the 
freeze out temperature lowered to 100 MeV 



Glauber, n/s = 0.0001, b = 7 f m 
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FIG. 2. (Color online) Pion V4,/(v2) 2 for ideal hydrodynamics 
calculations with varying freeze out temperatures. 
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FIG. 3. (Color online) V4,/(v2) 2 as a function of pt for identi- 
fied particles from ideal hydrodynamics with Glauber initial 
conditions and 140 MeV freeze out temperature. 



of Rcf. 2]). Indeed, CGC initial conditions, with a larger 
initial eccentricity, lead to a smaller value than Glauber 
while more peripheral collisions have a larger value than 
more central collisions, despite also having a larger ec- 
centricity. 

The effect of raising the freeze out temperature to a 
more realistic value can be seen in Fig. [2] The choice 
of freeze out temperature makes a significant difference 
to the shape of the ^/(i^) 2 curve. This sensitivity to 
Tf is perhaps not a complete surprise, since it is known 
from transport calculations that adjusting the treatment 
of freeze out does affect v 4 [H- For T f = 140 MeV, 
Vi/iy^) 2 is very fiat as a function of pt (reminiscent of 
data), which is not the case when Tf is significantly larger 
or smaller. It is interesting that this freeze out temper- 
ature is the same as gives the best fit of charged hadron 
multiplicity, (pt), and i>2 3 in viscous hydrodynamics. 
Note that in Fig. [3J viscosity has been neglected and so 



it corresponds to a vi that is larger than seen experimen- 
tally. Nevertheless, recall that the initial eccentricity — 
and by extension the magnitude of i>2 — does not signifi- 
cantly effect v^/(v2) 2 , and so this coincidence may indeed 
be significant. 

Figure |3] shows the results from ideal hydrodynamics 
for identified particles when the freeze out temperature 
is raised to 140 MeV. Data show approximately the same 
value for identified particles in the measured momentum 
range [It], Gil ■ This is confirmed in our calculation for 
Pt ^ 1-1.5 GeV, below which it begins to depend on 
the parameters of the simulation as well as increasing 
numerical error (especially for the heavier particles for 
which V2 ~ 0). 

Finally, Fig. 2] shows the results of viscous hydrody- 
namic simulations for several impact parameters, with 
ideal hydrodynamics results included for comparison. 
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FIG. 4. (Color online) Charged hadron t>4/(i>2) 2 as a function of pt from ideal and viscous hydrodynamics with Glauber and 
CGC initial conditions. Experimental results have been scaled down to account for the estimated effect of fluctuations [2p| ]. 
Filled and open squares are for charged hadrons at 10-20% (with scale factor K — 1.56) and 40-50% (K — 1.38) centrality, 
respectively, from STAR [11]. Filled and open triangles are for charged hadrons at 15-20% (K = 1.47) and 45-50% (K = 1.38) 
centrality, respectively, from PHENIX [191 ] . Viscous results correspond to the bold lines in Table HI 



Also shown is experimental data from STAR [l5| and 
PHENIX [19j . Since eccentricity fluctuations increase 
V 4 /(v 2 ) 2 by a factor K = (e 4 )/(e 2 ) 2 0, we have es- 
timated K in each centrality bin using the PHOBOS 
Glauber Monte Carlo (23|, and scaled down the data by 
K. This scaling removes most of the (already small) cen- 
trality dependence for both STAR and PHENIX in this 
centrality range: the residual centrality dependence is 
much smaller than in our ideal hydrodynamics results. 
The larger values of ViKv^) 2 for STAR as compared to 
PHENIX can most likely be attributed to larger nonflow 
effects [H. 

It turns out that different choices for the form of 
the viscous correction Sf to the distribution function at 
freeze out can significantly change the shape of these 
curves (see Ref. [24j). Several aspects seem to be uni- 
versal for reasonable choices of Sf, however. In contrast 
to expectations [25] as well as the results at low freeze out 
temperature (not shown), viscosity does not universally 
increase the ration of i>4/(i>2) 2 - In fact, in most cases at 
this realistic freeze out temperature, the value was found 
to decrease compared to ideal hydrodynamics. Rather 
than making the agreement with data better, viscosity 
tends to make it worse. In particular, the correction in- 
duced by viscous effects depends strongly on centrality: 
viscous effects are larger for the smaller systems created 
in peripheral collisions. It should be noted that choos- 
ing a Sf with a weaker momentum dependence than the 
standard quadratic ansatz (used in Fig. |4] as well as all 
previous viscous hydrodynamics results) can flatten these 
curves as well as reduce the strong impact parameter de- 
pendence (even potentially reducing it beyond that of 
ideal hydrodynamics), and this could in principle provide 
information about the dynamics of the hadron resonance 
gas at freeze out [24j |, 



V. LHC PREDICTIONS 



Using what we have learned from RHIC, it is then pos- 
sible to predict what will be seen in heavy ion collisions 
at the LHC. Using the same values for tq, Tf, and r]/s 
that led to the best fit of RHIC data 0, T% was adjusted 
to reproduce the expected multiplicity at top-energy Pb- 
Pb collisions at the LHC, as were the appropriate nu- 
clear parameters (the same as for the LHC predictions 
in Ref. [2a |. to which we refer the reader for details. For 
reference T { is 420 MeV and 390 MeV for Glauber and 
CGC initial conditions, respectively.) 



The results are shown in Fig. [5j LHC results for ideal 
hydrodynamics are similar to RHIC results with a lower 
freeze out temperature (see Fig. [2]). Specifically, there 
is an increase at small pt and viscous corrections can be 
positive as well as negative. However, here the viscous 
corrections are much smaller than for realistic RHIC sim- 
ulations. This means that — unlike for RHIC — the choice 
for the viscous correction to the non-equilibrium distribu- 
tion function (Sf) is relatively unimportant, and a com- 
paratively robust prediction can be made, with just a 
correction factor expected from fluctuations put in by 
hand. For example, a measurement of v^/(v2) 2 in a cen- 
trality bin of 20-60% is expected to be larger by a factor 
of ^1.34. Although there is some dependence on initial 
conditions and impact parameter, the results are very 
similar to those obtained for RHIC, so it is reasonable to 
make the prediction that ^/(i^) 2 measured at the LHC 
will be similar to that measured at RHIC, except for a 
small increase with decreasing p t . 
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Glauber CGC 




Pt (GeV) p, (GeV) 

FIG. 5. (Color online) Charged hadron V4,/{v2) 2 for 5.5 TeV Pb-Pb collisions at the LHC. The same parameters were used as 
for the predictions in Ref. [2a |. The measured ratio of separately averaged quantities should be larger by a factor of ~1.25-1.4, 
depending on the centrality bin (K = 1.34 for a bin at 20-60% centrality). 



VI. CONCLUSION 



In summary, we presented calculations of va/{v2) 2 from 
ideal and viscous hydrodynamic simulations of heavy ion 
collisions at RHIC, as well as those planned at the LHC. 
We confirm previous results from low Tf ideal hydro- 
dynamics and transport calculations. This observable 
is found to be insensitive to initial eccentricity, but to 
depend strongly on the freeze-out temperature Tf. In 
ideal hydrodynamics, the dependence on p t is found to 
be flat— as in RHIC data— only if T f is near 140 MeV, 
which is precisely the value providing the best fit to pt 
spectra and i>2. In addition, ui/(i>2) 2 was found to be the 
same for pions, kaons, and protons — also in accord with 
data. There is an unexpected dependence on impact pa- 
rameter, however, that is not seen in data. Including vis- 
cosity using the standard freeze out scheme causes both 



the shape of the curve as a function of p t and the depen- 
dence on impact parameter to deviate from what is seen 
in data (although both aspects could potentially be fixed 
by using a more appropriate form for the viscous cor- 
rection to the equilibrium distribution function at freeze 
out— see Ref. dj). 

For the LHC, this uncertainty in the freeze out scheme 
is less important, and we find that ^4/(^2) 2 should be 
similar to the value at RHIC, with a slight increase at 
small transverse momentum. 
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